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We investigate by means of ab initio electronic structure 
analysis and Quantum Monte Carlo calculations the sce- 
nario where longer-ranged magnetic interactions dominate 
over short-ranged interactions in the physical description of 
compounds. This question is discussed, in particular, for the 
case of CaCuGe2 06 which shows a spin-singlet behavior in- 
duced by third nearest neighbor copper pairs. 

PACS numbers: 75.30.Gw, 75.10.Jm, 78.30.-j 



There are a few low-dimensional materials which show, 
contrary to initial intuitive arguments, longer-ranged 
dominated exchange couplings. Examples include sys- 
tems like the most discussed 1/5-depleted Heisenberg 
system CaY^Oga or the alternating chain compound 
(VO)2P207D. CaV40g which was originally viewed as an 
array of weakly coupled square plaquettes of V 4+ ions, 
each in a singlet in the ground state, is found to be a sys- 
tem with strongest coupling between next nearest neigh- 
bor vanadiumaj as revealed by a recent electronic struct 
ture studyo. In (VO)2P2C>7 the strongest exchange pathEl 
is found to be that between two V 4+ ions through two 
phosphate groups PO4 and not between nearest .neighbor 
vanadium ions V-O-V as was initially thoughta. From 
these examples, we learn that the subtle competition be- 
tween the various interactions can lead to surprises and 
a microscopic study is needed. The detailed investiga- 
tion of these systems often unveils interesting properties. 
One such a system is CaCuGe206,i-a system related to 
the spin-Peierls compound CuGeOatl 

Susceptibility measurements!!] on CaCuGe2 0g show 
the existence of a spin-singlet ground state with an en- 
ergy gap of 6 meV. But, contrary to CuGeO-3 where, the 
spin-gap opens at the spin-Peierls phase transition^, in 
CaCuGe2 0g the spin gap is intrinsic and there is no phase 
transition between 4.2K and 300K. Also, inelastic neu- 
tron scattering (INS) measurements were carried out on 
CaCuGe2 0g powder and the existence of a finite spin-gap 
was confirmedcl. 

A question which neither INS nor susceptibility mea- 
surements could clarify was which Cu-pairs form the sin- 
glets in this material. This issue is of great importance 
since it determines the magnetic properties of the sys- 
tem. The structure of this material shows an obvious 
zig-zag ID chain of spin=l/2 Cu 2+ ions along the c-axis 
(see Fig. |l|). Nevertheless, the magnetization and suscep- 
tibility data of Sasago et aln are in disagreement with a 
spin=l/2 Heisenberg chain model. These authors pro- 



posed that the spin-gap behavior of the system should 
be due to longer-distance-formed Cu-pairs and consid- 
ered possible scenarios depending on which Cu-pairs were 
taken into account. Also the analysis of INS powder 
data led Zheludev et aln to conclude that this compound 
should be described as a weakly interacting dimer system 
in spite of the fact that the material shows a pronounced 
ID arrangement of magnetic ions. Which Cu ions form 
the pairs, however, couldn't be strictly defined from the 
analysis of the above experiments. 

It is the purpose of this letter to investigate in detail 
the singlet nature of this system and why it shows such 
a behavior. In order to do that we have performed (i) 
ab initio calculations which give us information on the 
chemical bonding and the electronic structure and (ii) 
Quantum Monte Carlo computations which are used in 
order to calculate thermodynamical quantities like the 
susceptibility and magnetization in this material. 

Crystal structure.- CaCuGe206 crystallizes in the 
monoclinic space-group P2i/c with lattice parameters a 
= 10.198 A,b = 9.209 A, c= 5.213 A, (3 = 105.73°, ar 
it contains four formula units per primitive unit celll 
The copper ions are all equivalent and form Jahn- Teller 
distorted octahedra with the surrounding oxygens and 
they build zig-zag chains along the c-direction (see Fig. 
H(a)) by sharing an edge of the octahedra. Intercalated 
between these chains are the Ge and Ca ions. The dis- 
tances between two nearest neighbor (NN) Cu within the 
same chain is 3.072 A and the Cu-O-Cu angles in these 
chains are 92° for the Cu-OIA-Cu angle, and 98° for the 
Cu-OIB-Cu angle (see Fig. || and Ref£3 for the oxygen 
notation). The Cu-Cu distance within the 2nd NN pair 
which are also located in the same zig-zag chain is 5.213 
A. The 3rd and 4th NN are formed by Cu atoms be- 
longing to neighboring be planes (see Fig. [j]). The dis- 
tance between two neighboring be planes alternate be- 
tween d = 4.4621 and d* = 5.354A giving rise to 3rd NN 
and 4th NN distances of 5.549 A and 6.213 A respectively. 

Band-structure.- We have performed an ab initio 
study based on the density-functional theory j-ODFT) in 
the generalized gradient approximation (GGA)li 4 ] in order 
to derive the electronic properties of CaCuGe2 06- Cal- 
culations have been performed within the frame-work of 
both the full-potential linearized augmented plane wave 
(LAPW) method based on the WIEN92L3 code and the 
linearized muffin tin orbital (LMTO)Ej method based 
on the Stuttgart TBLMTO-47 code. The results of the 
band-structure calculations show four narrow bands close 
to the Fermi level with a bandwidth of ~ 0.5 eV, having 
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predominant Cu d x 2_ y 2 character (in the local frame) ad- 
mixed with oxygen contributions. These bands are half- 
filled and strong correlation effects should explain the 
insulating groundstate in this compound. These bands 
are separated by an energy gap of about 0.3 eV from 
the valence-band set and a gap of about 2 eV from the 
conduction bands (see Fig. g). 

Effective model.- We have performed LMTO-based 
downfolding and tight-binding analysis to derive the 
(single-particle) tranfer integrals of an effective model for 
this system. 

In recent years a new form of the LMTO method has 
been proposed and implementedt^l which allows to de- 
rive few-orbital effective Hamiltonians by keeping only 
the relevant degrees of freedom and integrating out the 
irrelevant ones. This procedure amounts to putting the 
inactive orbitals by downfolding in the tails of the ac- 
tive orbitals kept in the basis. As a result, this takes 
into account the proper renormalization in the effective 
interactions between the active orbitals. The effective in- 
teractions obtained are therefore unique and given by the 
intervening interaction paths. Application of the down- 
folding procedure to CaCuGe2C>6 to derive one orbital 
(namely d x 2_ y 2 per Cu atom) effective model leads to the 
unambiguous result that the singlet pairs are formed by 
3NN Cu-pairs. The next important interactions are the 
INN providing a picture of two-dimensional interacting 
dimers in agreement with susceptibility data and INS. It 
is to be noted here that a simple-minded tight-binding fit- 
ting method with one-band model is unable to distinguish 
the relative importance between 3NN and 4NN pairs. Ei- 
ther of the choice of dimer pairs would lead to similar 
band dispersions though not to the same magnetic behav- 
ior. This dominant effect of the 3NN Cu-Cu interaction 
over the 4NN Cu-Cu interaction is originated from the 
specific orientation of the intervening oxygen atoms be- 
tween the Cu pairs and their relative hybridization with 
Cu d x 2_ y 2 orbital guided by the magnitude of Cu-O-Cu 
bond angles and the bond lengths. This is supported by 
the 3D electron density plot shown in Fig. || near the 
Fermi level in a plane perpendicular to the c-axis. As 
it can be seen, the exchange paths connecting two 3NN 
coppers are not equivalent to the exchange paths con- 
necting two 4NN coppers. For the considered 3NN Cu- 
Cu nait,the intervening oxygens (predominantly 02A and 
03AErt9) provide a well-defined interaction path while for 
the 4NN Cu-Cu pair the contribution of the intervening 
oxygens (predominantly 02B and 03B) is much smaller 
leading to almost disconnected Cu 2+ ions. This fact is 
also supported by the partial density of states (DOS) 
analysis at the Fermi level of the oxygens' contribution 
where we observe that the contribution coming from 02B 
and 03B orbitals is almost negligible compared to that 
of 02A and 03A. 

The downfolded-TJi bands arc represented by a Hamil- 
tonian of the form: 



H TB = J2u j (4c j + h.c.) (1) 

where the tij are hopping matrix elements between Cu 
neighbors i and j. We have considered hopping matrix 
elements up to fourth nearest neighbor (see Fig. |l|) t±, i2, 
t 3 and ti (the subindex denotes INN, 2NN, 3NN and 4NN 
respectively). The values derived for this effective model 
were found to be ii=0.068eV, t 2 =0.008eV, t 3 =0.088eV 
and i4=0.004eV. In the inset of Fig. (|J) we show the com- 
parison of the DFT-bands to the downfolded-TB bands 
considering t±, ti, £3 and £4 hoppings only. The band- 
splitting along the path DZ is reproduced by considering 
longer-ranged hopping matrix elements. 

An estimate of the exchange integral related to the 
most dominating 3NN interaction parameter £3 can be 
obtained by using the relation J 3 ~ 4i 2 /£7 where U is 
the effective on-site Coulomb repulsion. A value of U 
~ 4.2eV was proposed for CuGeC-3 bit-jnapping experi- 
mental data onto one-band descriptiontj. Assuming that 
this value is similar for CaCuGe206, we get a value of 
J 3 ~ 86K. 

Susceptibility and Magnetization.- The descrip- 
tion for CaCuGe20g obtained from the previous first 
principle calculation is that of a system of interacting 
dimers. In order to check this result we have ana- 
lyzed both the susceptibility as well as the magnetiza- 
tion behavior of this material by considering the following 
model: 

H e ff = ^3 ^ SiSj + J\ ^ S i$j > ( 2 ) 

<ij> 

where J3 and J\ are the exchange integrals between 3NN 
and INN as illustrated in Fig. (||). The analysis of the 
model has been done by Quantum-Monte-Carlo (QMC) 
simulations (stochastic series expansionEZMla) on 20 x 40 
lattices. We found an optimal value for J3 = 67 K, 
very close to the value 68 K proposed by Sasago et alB. 
We find clear evidence for a weak ferromagnetic inter- 
dimer coupling J\ < 0, as shown in Fig. (0). The op- 
timal value is J\ = — O.2J3. The g-factor was deter- 
mined by fitting the QMC-results for the susceptibil- 
ity, — {{S z — (S z )) 2 ) with the experimental sup. 
ceptibility x ( m [emu/mol]) at high temperatures viat2l 
X = 0.375(<7 2 /J)x^' 1 ^. The optimal g-factor was found 
to be g — 2.032, close to the value 2.07 proposed earlier!! 

Sasago et al. measured the magnetization curve for 
CaCuGe20g as a function of an applied magnetic field at 
T = 4.2 K. With the stochastic series expansion QMC- 
method it is possible to simulate quantum-spin models 
in an external fieldta, although these simulations become 
very hard near saturation. We present the results of this 
simulation in Fig. (||), for the same parameters we used 
for evaluating the susceptibility, together with the exper- 
imental results. It is very reassuring that J\ = — 0.2 J3 is 
again the optimal parameter set. 

The phase diagram of model Eq. (^|) is very interesting 
by itself. In the limit J3 = it consists of decoupled J\- 
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chains (see Fig. (0(a))) either antiferro- or ferromagnetic 
for Ji > or Ji < respectively. In both cases the 
excitation spectrum is gapless, while it shows a gap in 
the limit J\ = J3 ^ (dimer phase). We therefore 
expect two quantum critical points connecting these two 
regions. We have investigated the phase-diagram and 
found the critical coupling strength to be J\ ~ O.55J3 
and Ji ~ — O.9J3 respectively. The parameter range for 
CaCuGe206 is away from these two critical points. 

Discussion and Conclusions.- the analysis of the 
electronic structure of CaCuGe 2 06 by first principles 
calculations as well as the examination of susceptibility 
and magnetization data by the QMC method leads to a 
unique description of this material as a system of dimers 
formed by 3NN magnetic spin=l/2 Cu 2+ ions with ferro- 
magnetic (INN) interdimer couplings. The ferromagnetic 
nature of the INN coupling is more subtle-tp understand 
than in other edge-sharing Cu-0 systemsEj since due to 
the strong distorsion of the Cu octahedron, this cou- 
pling proceeds through two pathes with angles 92° and 
98°. These, two pathes can provide in principle compet- 
ing ferro-Eil and antiferromagnetic couplings respectively. 
The analysis of the susceptibility and magnetization be- 
havior shows nonetheless a clear effective ferromagnetic 
coupling. In this system, since the primary role is played 
by the 3NN Cu-pairs, the role of the possible frustrat- 
ing 2NN term to the INN interaction is secondary, un- 
like the case of the related spin-Peierls material CuGeC>3. 
The strength of the magnetic exchange integrals has been 
obtained from the analysis of the susceptibility and the 
magnetization data with very good agreement between 
these two measurements. DFT calculations overestimate 
the value of J\ , which is not surprisingc3, nonetheless the 
physical picture of this material is correctly provided by 
DFT + downfolding-TB analysis. 

In conclusion, by analizing the electronic and magnetic 
properties of CaCuGe20g we have been able to provide a 
physical picture of the system as well as a microscopic de- 
scription of the singlet nature in this material. The com- 
bination of ab initio techniques with many-body methods 
(QMC) can unambigously determine the microscopic be- 
havior of a system as has been shown for CaCuGe206 in 
particular. 

CaCuGe206 shows an interesting spin-gap behavior 
and we believe that experiments under the application 
of a magnetic field would be desirable to investigate the 
possible appearance of long-range order in the system. 
Moreover, the 2D interaction model provided by this sys- 
tem Eq. (j^) is interesting per se since it shows two quan- 
tum critical points as a function of the strength quocient 
J1/J3. 
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FIG. 1. (a) Projection on the be plane of the Cu zigzag 
chains in CaCuGe 2 6 , where the INN, 2NN and 3NN Cu 
pairs are indicated. The chains alternate between two neigh- 
boring be plane (denoted by black and grey balls) . (b) Projec- 
tion on the ab plane of the Cu sites, where the 3NN and 4NN 
Cu pairs are indicated. The linked Cu sites correspond to 
INN. d and d* denote the distances between two neighboring 
be planes. 
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FIG. 2. DFT results for CaCuGe 2 6 . The path is along 
r=(0,0,0), B=(-7T,0,0), D=(-7r,0,7r), Z=(0,0,tt), T, Y=(0,7r,0), 
A=(-7r,7r,0) E=(-7r,7r,7r). Also shown in rectangles is the Cu 
band character (in the local coordinate system) a,— d x 2_ y 2, 
b=d 3z 2_ 1 , c=d xy ,d xz ,d yz Shown in the inset are the four 
DFT-bands close to the Fermi level (left panel) and the cor- 
responding downfolded-TB bands considering up to 4-th NN 
hopping (right panel) along the same path. 
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FIG. 3. Electron density plot (right panel) for bands close 
to the Fermi level in a plane perpendicular to the c-axis. Note 
the dominant electron density for the 3NN path with respect 
to the 4NN path. The left panel shows the crystal structure in 
the corresponding plane. Big balls represent Cu atoms while 
small baAtei represent O atoms. The various non- equivalent 
oxygensulij are marked in the figure. 
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FIG. 4. Temperature dependence of the magnetic sus-i 
ceptibility for CaCuGe2 06- Shown are the experimental!] 
data (stars) and the QMC results for the model Eq. (^) 
for ,h = 67K and Ji/J z = 0.2 (up-triangles, g = 2.076), 
J1/J3 = 0.0 (open circles, g = 2.032), J1/J3 = -0.2 (full line, 
g = 2.032) and J1/J3 = -0.4 (down-triangles, g = 2.011). 
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FIG. 5. Magnetization fot, CaCuGe 2 6 at T = 4.2 K. 
Shown are the experimental!] data (stars) and the QMC 
results for the model Eq. (^) for J1/J3 = 0.0 (open cir- 
cles), J1/J3 = —0.2 (full squares) and J1/J3 = —0.4 
(down-triangles) with the same parameters as in Fig. (Gj). 
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